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1 . INTRODUCTION 


This  interim  report  provides  a summary  of  the  technical  activities 
pursued  under  contract  F30602-76-C-0097  with  RADC  during  January-December 
1976.  Sections  2 and  3 contain  a brief  description  of  the  work  on  Bayesian 
software  correction  limit  policies  and  a software  reliability  model  with 
imperfect  debugging,  respectively.  Section  4 lists  the  current  work  and 
plans  for  the  next  period. 
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2.  BAYESIAN  SOFTWARE  CORRECTION  LIMIT  POLICIES 


This  study  deals  with  the  problem  of  determining  an  optimum  correction 
limit  policy  for  a large  software  system  which  is  subject  to  random  occur- 
rences of  errors.  When  an  error  occurs,  a corrective  action  is  undertaken 
to  remove  it.  Such  an  action  can  be  scheduled  at  two  levels,  which  we  call 
Phase  I and  Phase  II.  By  Phase  I we  mean  that  the  corrective  action  will 
be  undertaken  by  the  programmer  while  Phase  II  action  is  undertaken  by  a 
system  analyst  or  system  designer.  First,  Phase  I corrective  action  is 
scheduled  for  a specified  time  T.  If  the  error  is  not  corrected  in  this 
time,  it  is  referred  to  Phase  II.  Figure  2.1  shows  the  sequence  of  corrective 
actions  in  an  operational  phase.  The  objective  is  to  determine  the 
optimum  value  T*  of  T which  minimizes  the  long  run  average  cost.  Two 
models  are  developed  for  this  purpose.  In  the  first  model  (Model  1) 
we  assume  that  the  cost  of  observations  of  error  occurrence  and  correction 
time,  prior  to  the  implementation  of  the  optimum  policy,  is  negligible. 

The  second  model  (Model  2)  incorporates  the  cost  of  observations. 

In  the  following  sections  we  provide  a brief  description  of  the 
various  aspects  of  this  problem.  A detailed  discussion  will  be  given 
in  a forthcoming  technical  report. 

2. 1 Model  Development;  Negligible  Cost  of  Observations  (Model  I) 

2.1.1  Assumptions 

In  the  development  of  the  model,  the  following  assumptions  are  made. 

(i)  The  error  occurrence  time  in  a software  system  has  an 
exponential  distribution  with  an  unknown  mean  X. 

(ii)  The  error  correction  time  at  Phase  1 is  exponential 
with  an  unknown  mean.*). 

(iii)  The  Phase  H error  correction  time  has  a general  distribution 
with  a known  mean  p2> 

(iv)  Reasonable  prior  distributions  can  be  chosen  for  X and 


2 


time 


Figure  2.1.  Sequence  of  Corrective  Actions 
in  Operational  Phase. 
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2.1.2  Predictive  Distributions 


Finally,  let  x = (x. , ....  x ) and  y = (y. , ....  y ) be  the  observed 
i n l n 

values  of  n error  occurrence  times  and  n error  correction  times,  re- 
spectively. 

Now  we  obtain  expressions  for  the  predictive  distributions  of  Y 
and  X and  the  Bayesian  estimates  ?n+^  and  x^+^  which  we  will  use  to 
obtain  the  cost  function. 

lor  given  observations  the  likelihood  function  of  y^  is 

n » 

- £ 

KUjJjr)  * -e  i_1  , 

and  the  posterior  distribution  of  y^  is 


4 


n a,+n 


P(Uj \y) 


+ £*!>  1 -(Oj+n+1)  *{(e1+1^y1)/(lJl)1 


roij+ro 


Then,  the  predictive  distribution  of  error  correction  time  at 
Phase  I is 


g(y|i.)  = (~ 


n + a. 


-)  ( 1 + 


-(n+a^+1) 


Z y + 6 
i=l 


E y + 6 
i=l 


and  the  cumulative  predictive  distribution  to  some  specified  time  t is 


G(t|^)  = 


g(y li)dt  = i - (i  + 


-) 


-(n+a1) 


E y + B 

i=l 


We  define  the  predictive  Phase  I error  correction  rate  as 

g(t  lx) 

r(t|£)  = — 


so  that 


G(t|y) 
n + a 


1 


r (t  |jr)  = (— ) (1  + 


-) 


E y + B 
i=l 


E y + Bx 

i=l 


where 

G(t)  = 1 - G (t ) 

From  the  above  results  the  Bayesian  estimate  of  the  (n+l)st  error 
correction  time  for  given  £ is 

- .kill 

^n+1  + n - 1 

Proceeding  similarly,  the  Bayesian  estimate  of  the  time  to  (n+l)st  error 
occurrence,  for  given  x,  is 

n 

E x f B2 

i=l  * 

it  n . 1 m m 

n+1  «2  t n - 1 
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2.1.3  Coat  Function 


Let  c^(c2)  be  tbe  cost  per  unit  time  of  error  correction  In  Phase  I 
(Phase  II)  and  the  costs  be  linear  functions  of  time.  If  we  consider 
one  cycle  to  be  the  time  from  the  beginning  of  (n+l)st  operation  to  the 
beginning  of  (n+2)nd  operation,  then  the  expected  cost  In  one  cycle  Is 

T_ 

E(C)  - c1j  G(tlx)<it  + c2  u2  G (T  |x> » 


where  T denotes  the  scheduled  correction  limit  time  in  Phase  I. 

The  expected  length  of  one  cycle  is 

T 

E(T)  - £n+1  + / G(t|Z)dt  + p2  G(T|x)  . 

0 

and  hence  the  long  run  exp  ;tc  cost  per  unit  time  Is 


cm  = 

‘ E(T) 


or 


C(T)  - 


ci  /G(t|i)dt  + c p2  G(T|^) 

0 

*n+l  + f5(tlz)dt  + H2  G(T|Z) 


2.1.4  Optimum  Policy 

From  the  above  expressions,  we  note  that 

c2  W2 


C(0) 


Xn+1  + W2 


and 


C(.)  - ^ 


Xn+1  + yn+l 
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where  y , , is  the  Bayesian  estimate  of  y for  given  data  y.  Also,  note 
n+l 

that  T = 0 means  that  the  errors  are  corrected  only  at  Phase  II  while 
T - oo  means  that  they  are  corrected  at  Phase  I. 

To  obtain  an  optimum  T*  which  minimizes  the  long  run  average  cost  per 
unit  time,  C(T),  we  need  the  following  theorems  and  corollary.  They  are 
given  here  without  proofs. 

Theorem  2.1: 

Assume  cj<c2’  Then  there  exists  a finite  and  unique  T*  which  satisfies 


r (T |jr)  {c2  xn+1  + (c^c^  J G(t|^)dt} 


+ (c2-c1)  G(T|x)  - 


c,  x ■ , 
1 n+l 


under  the  following  condition 


Cl(*n+1  * ^ ~ C2 

C2  Vl  W2 


Theorem  2.2: 

If  the  above  conditions  are  satisfied  then  there  also  exists  a finite 
and  unique  upper  bound  T(>T*)  such  that 

C1  X J.1 

_/Si  x _ 1 n+l  


y2{c2  *n+l  + (C2-Cl)yl} 


This  upper  bound  can  be  used  to  obtain  an  Initial  value  for  solving 
the  nonlinear  equations  in  T*. 

Corollary  2.3 

If  there  exists  an  optimum  T*  then  the  associated  cost  function  is  given 


C(T*)  - 


Cl"C2  M2  r^T*l^ 


1 - P2  r(T*|i) 
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2.1.5  Numerical  Example 


We  use  simulated  data  in  this  example  to  illustrate  the  calculation 
and  nature  of  various  quantities  in  the  determination  of  T*. 

Let 

cL  - 8000  c2  - 9000 

- 0 B - 0 

a2  “ 0 B2  “ 0 

M2  - 0.7 

The  simulated  data  (x  ,y  ) are  given  in  Table  2.1.  Suppose  n - 10  data 

n n 

points  are  available.  The  Bayesian  estimates  of  failure  time  and  correction 
time  are  x^  - 59.60  and  y^  - 0.78,  respectively.  Such  values  for  various 
n are  given  in  Table  2.2  . For  the  case  n - 10  we  see  that  the  optimum 
correction  limit  time  is  T*  ■ 0.90  hours  and  the  corresponding  minimum 
cost  rate  is  C(T*)  - 99.44  dollars/hour. 

Thus,  for  this  set  of  data,  we  will  schedule  corrective  action  in 
Phase  I for  0.90  hours  and  if  it  cannot  be  completed  in  this  time,  the 
software  system  will  be  referred  to  the  system  programmer  for  corrective 
action. 
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Table  2.1 


Simulated  Values  of  x and  y 

n 


n 

x (Hrs.) 
n 

yn  (Hrs.) 

n 

x (Hrs.) 
n 

yn(Hrs 

i 

61.34 

1.90 

11 

53.44 

1.03 

2 

27.84 

1.08 

12 

2.87 

0.95 

3 X 

154.30 

0.85 

13 

31.27 

0.60 

4 

14.58 

0.26 

14 

97.06 

0.02 

5 

10.86 

0.01 

15 

78.17 

1.49 

6 

35.35 

0.31 

16 

124.52 

0.52 

7 

140.13 

0.38 

17 

0.49 

0.36 

8 

36.47 

1.50 

18 

12.33 

0.08 

9 

8.74 

0.43 

19 

85.44 

3.51 

10 

46.79 

0.27 

20 

23.59 

0.10 
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Table  2.2 


Calculation 

for  the  Optimum 

Correct  Ion 

Time  Policy 

n 

Vi(hr>) 

Vi(hr'> 

T*  (hr.) 

C (T*  ) 

2 

89.17 

2.98 

0 

70.10 

3 

121.74 

1.92 

0 

51.45 

4 

86.02 

1.36 

0 

72.65 

5 

67.23 

1.02 

0 

92.74 

6 

60.85 

0.88 

0.32 

101.05 

7 

74.07 

0.80 

0.73  " 

80.11 

8 

68.70 

0.90 

0.02 

90.78 

9 * 

61.20 

0.84 

0.38 

100.60 

10 

59.60 

0.78 

0.90 

99.44 

11 

58.98 

0.80 

0.66 

102.87 

12 

53.88 

0.82 

0.50 

113.82 

13 

52.00 

0.80 

0.68 

116.85 

14 

55.46 

0.74 

1.45 

103.80 

15 

57.09 

0.79 

0.75 

106.51 

16 

61.58 

0.77 

1.02 

97.47 

17 

57.76 

0.75 

1.45 

101.22 

18 

55.09 

0.71 

2.16 

101.14 

19 

56.78 

0.86 

0 

109.61 

20 

55.03 

0.82 

0.13 

112.97 
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2.2.  MODEL  DEVELOPMENT  INCLUDING  SAMPLING  COST  (MODEL  2 


In  this  section  we  develop  the  cost  model  for  the  pptiraum  correction 
limit  policy  by  incorporating  sampling  cost.  It  1b  assumed  that,  in 
addition  to  the  assumptions  of  Section  2.1.1,  the  sampling  cost  is  a 


linear  function  of  sample  size. 

Let  c be  the  sampling  cost  at  each  stage,  and  (T^)  be  the  expected 

cost  per  unit  time  until  the  completion  of  (n+l)st  corrective  action  under 

scheduling  time  T . If  we  decide  to  take  another  observation  ((n+l)st), 

then  C , (T  . , ) is  the  cost  rate  function  until  the  completion  of  (n+2)nd 
n+1  n+1 

corrective  action  under  scheduling  time  Tn+^. 

2.2.1  Cost  Function 

Following  a procedure  similar  to  that  of  Sections  2.1,2  and  2.1.3, 


the  expected  cost  per  unit  time  at  the  end  of  the  (n+l)st  corrective  action  is 


nc  + c / G(t)2>dt  + c m2  G(T  lz> 
l 0 *■  n 

Cn(Tn)  = — T 

1 Xi  + Xn+1  + 1 yi  + r G(t |X)dt  + U2  G(T ll) 

i-1  1 i-1  1 0 


If  we  decide  to  take  the  next,  (n+l)st,  observation,  then  the  cost  rate 
function  to  the  end  of  (n+2)nd  corrective  action  is  similarly  obtained  as 


(n+l)c+c1  / G(t|;jr)dt  + c,  W2  g(t|x> 

_0 

"*n " ii  ^n+1 

Xi  + 2 Vl+  J/i  + Vl  + / 5(t|X)dt+  M2  G(Tn+1|Z) 
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2.2.2  Optimum  Policy 

Our  objective  Is  to  determine  the  optimum  sample  size  n*  and  the 
optimum  correction  limit  T^*  such  that 

C (T  *)  < C . . (T*  . ) 
n n n+1  n+1 

Given  that  n observations  have  been  taken,  the  following  steps  summarize 
the  procedure  of  determining  these  quantities: 

(i)  Calculate  Cn(T*)  and  Cn+1(T*+1) 

(ii)  If  C (T*)  *■  C (T*  ),  then  stop  taking  observations  and 

n n n+i  n+1 

employ  n*  and  T*  as  the  optimum  policy. 

(lii)  If  C (T*)  > C ..(T*),  take  (n+l)st  observation,  i.e.  let 
n n n+l  n 

n - n+1  and  go  to  step  (i). 

The  following  theorems  are  useful  in  determining  cn(T*)  and  Cn+^(T*+^). 
These  are  given  here  without  proofs. 


Theorem  2.4 


Suppose  c^  < c^  and  A * 0.  Then  there  exists  a unique  and  finite  T* 


satisfying 


where 


r<Tn  li>  (A  + (Cj-Cj)  J G(t  |jr)dt}  + ( Cj-Cj ) G(T^  |jr)  - B 


A - c0  { I x + x . + Z y } - nc 

2 i-i  1 n+1  i-i  1 

1 r n . n -i 

B ■ — c,{  E x.  + x , + I y . } - nc 

U2  L 1 i-i  1 n+1  i-1  1 -I 


Also,  the  associated  cost  rate  function  is  given  by 


c - c2  p2  r (T*) 
“ 1 - »2  r (T*) 
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Theorem  2 . 5 


If  the  conditions  of  Theorem  2.4  hold,  then  there  exists  a finite  and 
_ * 

upper  bound  T (>  T ) such  that 
n n 


r(TnIZ) 


A +'  (c2~°1)9n+i 


Theorem  2. 6 


If  both  T*  and  T exist,  then  the  following  relationship  holds: 
n n+1 

VV  (>, 
r<T")  ,<)  '‘Vl’ 

* it 

<«>  T > T . . 
n (<)  n+1 


2.2.3  Numerical  Example 

In  this  example  we  use  simulated  data  to  illustrate  the  determination 
* * 

of  n and  T . Simulated  values  of  x and  y for  various  n are  given  in 
n n n 

Table  2.3.  Let 


c1  - 8000, 


9000, 


40 


&1-! 

e2  - o 


P2  - 0.7 


Then  the  calculated  x^,  Tq,  C^),  Tr+1  and  Cn+1tTn+1)  are 

computed  from  the  above  expressions  and  are  given  in  Table  2.4,  From  this 
table  we  see  that  for  n - 11,  (T*^)  - 21,74  and  C12^T12^  “ 23,63  80  that 
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Table  2.J 


Simulated  Values  of  x and  y 

n n 


n xn (hr ) yn(hr)  n *n(hr)  yn(*»r) 


0.09. 


* * * * 
C,,(T,,)  < C „ (T .) . Therefore,  the  optimum  policy  is  n =11  and  T = 
ii  ii  iz  iz  n 

2 • 3 Concluding  Remarks 

In  the  previous  sections  we  have  given  an  overview  of  the  key  results 
for  obtaining  correction  limit  policies  for  two  models.  For  the  first  model 
the  sampling  costs  are  assumed  to  be  negligible  while  for  the  second,  the 
sampling  cost  is  taken  to  be  a linear  function  of  sample  size.  The  main 
theorems  have  been  given  without  proofs  and  the  policy  determination  was 
illustrated  via  two  numerical  examples. 

A detailed  technical  report  on  this  task  is  being  completed  for  sub- 
mission to  RADC. 
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3.  A SOFTWARE  RELIABILITY  MODEL  WITH  IMPERFECT  DEBUGGING 


The  purpose  of  this  study  is  to  develop  software  reliability  models 
for  the  case  when  errors  are  not  corrected  with  certainty.  In  other  words, 
we  are  interested  in  developing  and  studying  software  reliability  models 
for  the  case  when  the  programmer  sometimes  fails  to  correct  a detected  error. 

A description  of  the  key  results  is  given  in  the  following  subsections. 

3. 1 Assumptions 

The  following  assumptions  are  made  for  developing  the  model. 

(i)  Errors  in  the  software  are  independent  of  each  other  and  have 
a constant  occurrence  rate  A. 

(ii)  The  probability  of  2 or  more  errors  occurring  at  the  same  time 
is  negligible.  Then,  if  there  are  n errors  in  the  software  at 
present,  the  distribution  of  time  to  the  next  failure  is 
f (t)  ” nA-e 

(iii)  The  correction  time  in  the  model  is  assumed  to  be  zero.  This 

assumption  may  not  hold  in  some  situations.  We  plan  to  develop 
another  model  which  will  include  error  correction  time. 

3.2  Model  Development 

Let  p be  the  probability  of  successfully  correcting  an  error, 

q be  the  probability  that  the  error  is  not  corrected,  q - 1-p, 

X(t)  be  a random  variable  denoting  the  syntem  state, 

N be  the  initial  number  of  errors. 

Let  and  PN  ^ represent  the  transition  probabilities  from  state  N to 

states  N-l  and  N,  respectively.  Then  the  cumulative  distribution  function  (cdf) 

of  the  time  to  next  error  at  state  N is 
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r i “NAt 

FN(t)  - 1 - e 

FN,N-1  ’ P 
PN,N  “ q 

In  general,  the  cdf  of  one  step  transition  for  the  underlying  Markov 
Process  is: 


where 


?1  (t ) - 1 - e 


Qij(t)  = PiJFi(t)* 
-lAt 


and 


ij 


OJ 


(P 

q 

f l 
0 


if  J = i-1 , 1 < i < N 

if  J * i 

otherwise 

if  j = 0 
otherwise 


This  can  be  represented  as: 


(Pij) 


N-l 

N 


1 2 
0 0 


0 

q 


N-l  N 


p q 

0 p 


The  above  expressions  constitute  the  basic  model  for  the  error  process 
with  parameters  N,  p and  X.  For  these  expressions  we  obtain  various  quantities 
of  Interest  in  the  following  subsections. 
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3 . 3 Distribution  of  Time  to  a Completely  Debugged  Software 


Let 

G (t)  = cdf  of  first  passage  time  from  i to  0 by  time  t, 
i f u 

Then,  by  conditioning  on  the  next  one  step,  we  get  the  renewal  equa- 
tions which  we  solve  by  Laplace-Stieltjes  transforms.  We  get  the  cdf  and 
pdf  of  first  passage  time  as 

CK.0(t>  * j.  V1  ' E"‘PXt> 

j=l 

and 

rn  o(t)  = l c]  JpA  e_:lpAt. 

j=l  J 

N 1-1 

where  (-1)J 

The  cdf's  of  the  first  passage  times  for  several  values  of  P,  N = 10, 
and  A = 0.02  are  shown  in  Figure  3.1,  As  expected,  the  cdf  for  a larger  p 
dominates  that  for  a smaller  p. 

3.4  Distribution  of  Time  to  a Specified  Number  of  Remaining  Errors 

In  many  situations  we  are  not  interested  in  a completely  debugged 
program  because  of  the  cost  involved.  We  may  be  willing  to  tolerate  a 
certain  number  of  remaining  errors,  say  np,  which  will  ensure  some  desired 
reliability.  The  distribution  of  time  to  nQ  is  then  of  interest. 
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CDF 


FIRST  PASSAGE  TIME 


Figure  3.1.  CDF  of  Time  to  Completely  Debugged  Software. 
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Usi  ,g  an  approach  similar  to  that  of  Section  3.3,  we  get  the  cdf 


and  pdf  as  follows. 


N-n„ 


G = l B4  {1  - e~(n0+j)pAt} 

N,n0  i-1  J,n0 


and 


N-n, 


8N  n (t)  = £ B1  n (n0+j)pAe  ° 

'0  1=1  0 


(nn+j)pXt 


where 


R N!  j-1  j 

J.nQ  nQ!  j ! (N-nQ-j)!  nQ  + j 


Plots  of  pdf  and  cdf  for  n^O,  1 , . . , 9 , A = 0.02,  N = 10  and  P = 0.9  are 
given  in  Figures  3.2  and  3.3,  respectively.  These  plots  are  self  explanatory. 


3 . 5 State  Occupancy  Probabilities 

In  this  section  we  develop  expressions  for  computing  the  estimated 
n : r of  remaining  errors  after  a specified  time  period, 

et 

PN  (t)  = P(X(t)  = nQ|X(0)  = N} 

’ 0 

Using  the  renewal  theoretic  approach  and  the  relationships! between 

P„  (t)  and  G„  (t),  we  obtain  the  following  results; 

N,nQ  N,n, 


0 

N-n, 


pN  n (t)  = l A.  { 

N,n0  k=l  *cinn 


-n^pt  -(nQ-k)pAt 


} for  0 < nQ  < N 


and 


PN,N(t)  " 6 


-NAt 
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PDF 
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TIME 

Figure  3.3,  CDF  o£  Time  to  a Specified  Number  (nQ) 
of  Remaining  Errors. 


where 


N! 

n ! k!  (N-n  -k)! 
n n 


Further,  an  estimator  of  the  number  of  errors  remaining  at  time  t is: 


X(t)  - l n0P  (t) 

n^<N  ° N*n0 

Plots  of  the  probability  distributions  for  nQ  - 0,1, 2,..., 10  when  N - 10, 
p - 0.9  and  A - 0.02  are  given  in  Figure  3.4. 

Now  suppose  we  want  an  estimate  of  the  number  of  errors  remaining  at 
time  t * 100.  We  have 
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P i 

N,v 

UOO) 

10 

0,152 

x 10"7 

9 

0.769 

x 10-6 

8 

0.175 

1 

O 

X 

7 

0.235 

x 10”3 

6 

0.208 

x 10-2 

5 

0.126 

1 

O 

X 

4 

0.530 

x 10_1 

3 

0.153 

2 

0.290 

1 

0.325 

Q 

0.164 

From  these  values  we  get 

X(100)  - l n P (100)  -1.7 
_ Sin  0 N,n 
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Figure  3.4.  Probability  Distribution  for  Given  Number  of  Remaining  Errors. 
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3 . 6 Expected  Number  of  Errors  Detected  by  Time  t 

Let  N(t)  be  a random  variable  denotine  the  total  number  of  errors 
detected  by  time  t.  Further,  let  M^(t)  be  the  expected  number  of  errors 
detected  by  time  t when  the  initial  number  of  errors  is  N,  i,e. 

^(t)  = E[N(t)  | X (0)  =N ] 

Using  Harkov  renewal  theory  and  Laplace-Stielt jes  transforms,  we 
obtain  the  following  results: 

, N tf-k+1 

MN(t)  = p I I B {1  - e-(k-l+j)pAt 

F k=l  j=l  J *k~1  } 

Plots  of  M^t)  for  N = 3, 6, 9,... 30,  X = .02,  and  p = 0.9  are  given 
in  Figure  3.5  and  are  self  explanatory. 

3.7  Concluding  Remarks 

In  the  previous  sections  we  described  the  key  aspects  of  an  imperfect 
debugging  model.  Expressions  for  various  quantitites  of  interest  were  de- 
veloped and  numerical  solutions  were  given  for  selected  values  of  unknown 
parameters  N,  X and  p. 

A detailed  technical  report  on  this  task  will  be  submitted  to  RADC 
in  Spring  1977. 

The  model  was  developed  on  the  assumption  of  zero  correction  time. 

For  the  follow-up  model,  we  will  include  the  tire  for  corrective  action. 

Also  to  be  investigated  is  the  problem  of  estimating  the  parameters  from 
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EXPECTED  NUMBER 


TIME 

Figure  3,5,  Mean  Number  of  Detected  Errors  versus  Time. 
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available  data.  In  a Bayesian  context,  we  plan  to  develop  models  when 


parameters  N,  X and  p have  prior  distributions. 

Detailed  technical  reports  on  this  task  will  be  submitted  to  RADC  as 
work  is  completed. 
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4.  CURRENT  WORK  AND  PLANS  FOR  THE  NEXT  PERIOD 


The  work  on  two  modelling  tasks  was  described  in  Sections  2 and  3. 

In  this  section  we  present  an  outline  of  the  work  being  pursued  and  planned 
for  completion  In  the  next  period. 

Further  Work  on  Software  Reliability  Modelling  with  Imperfect  Debugging 
The  following  aspects  of  this  problem  will  be  investigated: 

(i)  Development  of  a model  to  incorporate  error  correction  time. 

(ii)  Estimation  of  parameters  of  the  models  being  studied. 

(iii)  Development  and  study  of  a Bayesian  imperfect  debugging  model 
Expressions  for  various  quantitites  of  interest  will  be  developed 
and  numerically  studied. 

Bivariate  Software  Error  Model 

We  plan  to  develop  software  reliability  models  when  the  errors  in 
the  system  can  be  classified  as  serious  and  non-serious.  Various  quantities 
of  interest,  for  example  system  availability,  will  be  studied. 

Software  Error  Data  Analysis 

Some  preliminary  work  i6  being  pursued  on  the  development  of 
empirical  models  for  software  error  data.  We  plan  to  continue  this 
activity. 
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Demonatr.it Ion  and  Testing  Plans 

Software  reliability  demonstration  plana  will  be  developed  to 
provide  a tool  for  making  accept/reject  decisions  for  software  packages. 
Both  classical  and  Bayesian  plans  will  be  investigated. 

Sensitivity  Studies 

In  this  task  we  will  investigate  the  effects  of  changes  in  prior 
distributions  and/or  model  parameters  on  quantities  of  Interest.  The  study 
will  be  mostly  numerical  in  nature. 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


Unit 


lengtli 

mass 

time 

electric  current 
thermodynamic  temperature 
amount  of  substance 
luminous  intensity 

SUPPLEMENTARY  UNITS: 

plane  angle 
solid  angle 

DERIVED  UNITS: 
Acceleration 

activity  (of  a radioactive  source) 
angular  acceleration 
angular  velocity 

area 

denaity 

electric  capacitance 

electrical  conductance 

electric  Held  strength 

electric  inductance 

electric  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  of  electricity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velocity 

viscosity,  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


metre 

kilogram 

second 

ampere 

kelvin 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  per  kelvin 

newton 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

joule 

watt  per  steradian 

toule  per  kilogram-kelvin 

pascal 

watt  per  metre-kelvin 

metre  per  second 

pascal-second 

square  metre  per  second 

volt 

cubic  metre 
reciprocal  metre 
joule 


SI  PREFIXES: 


Multiplication  Factors 

1 000  000  000  000  * 10’1 
1 000  (XX)  000-  10’ 

1 000  000”  10* 

1 000  » 10’ 
too  - to1 
10  ” 10' 

0 1 ” 10-' 
001  - 10-» 
0001  ” 10-' 
0.000  001  ” 10- 4 
0.000  (KH)  001  ” 10-  4 

oooo  ooo  ooo  ooi  • io  " 
0 .000  000  000  QOO  001  ” IO'" 

ooooooooooooooooooi  - in-'4 

• To  be  avoided  where  possible 


SI  Symbol 


m 

kg 

s 

A 

K 

mol 

cd 


bormula 


rad 

sr 


F 


H 

V 

V 

I 

N 

Hz 

lx 

Im 

Wb 

T 

A 

W 

Pa 

C 

I 


Pa 


V 

I 


m/s 

(disintegration  ]li 

rad/s 

rad/s 

m 

kg/m 

A-s/V 

A/V 

V/m 

Vs/A 

W/A 

V/A 

W/A 

N-m 

j/K 

kg-m/s 

(cycle|/s 

I m/m 

cd/m 

cd-ar 

A/m 

V-s 

Wb/m 

j/» 

N/m 

A-s 

N-m 

W/sr 

I'kg-K 

N/m 

W/m-K 

m/s 

Pa-s 

m/s 

W/A 

m 

(wave|/m 

N-m 


Prefix  SI  Symbol 


tera 

gig* 

mega 

kilo 

hecto‘ 

delta* 

decl* 

cent  I* 

mllll 

micro 

nano 

pico 

femto 

alto 


T 

G 

M 

k 

h 

da 

d 

c 

m 

M 

n 

r 


l % 
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i MISSION  j 

$ 0/  * 

^ /?om£  v4tr  Development  Center  2 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  coimunlcatlons 
(C3)  activities,  and  in  the  C3  areas  of  Information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  coimunications , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology , 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 
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